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Summary. — Gamma-ray emission produced by interactions between cosmic rays 
(CRs) and interstellar gas traces the product of their densities throughout the Milky 
Way. The outer Galaxy is a privileged target of investigation to separate interstellar 
structures seen along the line of sight. Recent observations by the Fermi Large 
Area Telescope (LAT) shed light on open questions of the EGRET era about the 
distribution of CR densities and the census of the interstellar medium. 
The gradient of 7-ray emissivities measured in the outer Galaxy is significantly 
flatter than predictions from widely used CR propagation models given the rapid 
decline of putative CR sources beyond the solar circle. Large propagation volumes, 
with halo heights up to 20 kpc, or a flat CR source distribution are required to match 
the data. Other viable possibilities include non-uniform CR diffusion properties or 
more gas than accounted for by the radio/mm- wave data. 7-ray data constrain the 
evolution of the Xco ~ A'^(H2)/Wco ratio within a few kpc from the Sun. There 
is a signiflcant increase by a factor 2 from nearby clouds in the Gould Belt to the 
local spur. No further significant variations are measured from the local spur to 
the Perseus spiral arm. At the level of statistical accuracy provided by the LAT 
data, the most important source of uncertainty, often overlooked so far, is due to 
the optical depth correction applied to derive the column densities of Hi. Reliable 
determinations of the amount of atomic gas in the plane are key to better probe the 
properties of CRs in the Galaxy. 

PACS 95.85.Pw - Astronomical observations: 7-rays. 

PACS 98. 38. Am - ISM in the Milky Way: physical properties. 

PACS 98 . 70 . Sa - Cosmic rays. 
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1. Introduction 

Interstellar 7-ray emission is produced by interactions of high-energy cosmic rays 
(CRs) with the gas in the interstellar medium (ISM) and the soft interstellar radiation 
fields. Its observations carry information about CR properties in distant locations and 
it provides a tracer of the total interstellar gas densities to be compared with radio/mm- 
wave data: the 21-cm line of atomic hydrogen, Hi, and the 2.6-mm line of CO, used as 
a surrogate tracer of molecular mass. 

Open issues in the understanding of the interstellar 7-ray emission concern the iden- 
tification and spatial distribution of CR sources and the census of the ISM, notably the 
-'^co = ^(H2)/Wco conversion factor. The outer Galaxy is a privileged observational 
target since the Dopplcr shift of radio/mm- wave lines due to the Galactic rotation un- 
ambiguously locates the emitting clouds. There are two longitude windows with a steep 
velocity gradient leading to a good kinematic separation [1]. 

We reported analyses of recent measurements by the Large Area Telescope (LAT) 
on board the Fermi j-ray Space Telescope [2] for the second [3] and third [4] Galactic 
quadrants. The component separation based on likelihood fitting allowed us to extract 
the emissivities per H atom, ^hi, and per Wqo unit, qco, in several regions along the 
lines of sight as described in table 1. We refer the interested reader to the aforementioned 
papers for details about the analysis and we briefiy discuss here the implications of the 
results for the distribution of CRs in the Galaxy and the calibration of the Xqo ratio. 



Table I. - Regions seen toward the outer Galaxy and their approximate Galactocentric distances. 





second quadrant 
100° < I < 145° 


third quadrant 
210° < I < 250° 


Gould Belt 


8.5 - 8.8 kpc 




Local Arm 


8.8 - 10 kpc 


8.5 - 10 kpc 


interarm region 




10 - 12.5 kpc 


Perseus arm 


10 - 14 kpc 


12.5 - 16 


outer region 


> 14 kpc 


> 16 kpc 



2. The spatial distribution of CR densities 

Provided that the H i column densities are accurately measured from radio data, the 
emissivity qui, i-e. the 7-ray emission rate per hydrogen atom, directly relates to the 
average CR densities in each of the regions considered. The qui profile (integrated above 
200 MeV) as a function of Galactocentric distance R is shown in fig. 1 (left panel). 

CR densities appear to be uniform within 20% for comparable Galactocentric dis- 
tances in the two regions studied. There is a small decrease from nearby complexes 
in the Gould Belt(^) to the Perseus spiral arm, but no further significant gradient is 
observed beyond i? 2± 11 kpc. 

We assessed the uncertainties due to the optical depth correction applied to the Hi 
line intensities to derive A^(Hi). In fig. 1 (left panel) we have explored a range of possible 



(^) Our measurement of the local emissivity is compatible with expectations based on CR 
spectra as they are measured near the Earth [3, 5]. 
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Fig. 1. - Left: gni (integrated above 200 MeV) as a function of Galactocentric radius R for 
different Hi spin temperatures, as measured in tfie second (black) and third (red or gray) 
quadrants. Statistical errors are comparable with marker sizes. Right: Xco as a function of 
R. The error bars include the uncertainties due to the H i optical depth correction. The (blue) 
step function represents the model by Strong et al. [6], the (green) line that by Nakanishi and 
Sofue [7]. 



spin temperatures (Tg), showing that the uncertainties in the Hi mass, often overlooked 
in the past, dominate the errors in the emissivity determination. 

3. The calibration of molecular masses 

If subject to the same CR fluxes, a hydrogen molecule emits two times 7-rays as a 
hydrogen atom. We can therefore calibrate Xco as 9co/2 9hi- Our results, constraining 
^co over a few kpc from the solar system, are shown in fig. 1 (right panel). 

There is a significant increase by a factor ^ 2 from the nearby clouds of Cepheus and 
Cassiopeia in the Gould Belt to the local arm, whereas no further significant variations 
are observed up to i? ^ 14 kpc. The increase by one order of magnitude in the outer 
Galaxy proposed by Strong et al. [6] is not confirmed. The 7-ray estimates are also 
systematically lower than predictions by Nakanishi and Sofue [7] from virial masses. 

4. The cosmic-ray gradient problem 

The comparison between the measured 7-ray emissivities and expectations from CR 
propagation models has implications for the origin and propagation of cosmic rays in the 
Galaxy. We adopted the widely-used GALPROP propagation code [8]. A conventional 
model analogous to [6, 9] (solid line in fig. 2, left panel) predicts a gradient steeper than 
inferred from the LAT measurements. 

To alleviate the discrepancy one can increase the CR propagation halo or the radial 
scale of the CR source distribution. In fig. 2 (left panel) we varied the halo height(^), 
showing that halo heights from 10 kpc to 20 kpc are preferred if the CR source profile is 
that derived from pulsar and supernova remnant observations. In fig. 2 (right panel) we 



(^) And correspondingly the diffusion parameters to stay consistent with CR isotopic abun- 
dances measured at the Earth. 
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Galactocentric Radius (kpc) Galactocentric Radius (kpc) 

Fig. 2. - I as a function of galactic radius R: bow-ties represent our estimates for the second 
(black) and third (red or gray) quadrants. The curves give the predictions by some GALPROP 
models, on the left varying the height of the propagation halo from 1 kpc to 20 kpc (the solid 
line corresponds to 4 kpc, a value often assumed in past studies [9]), on the right setting the 
density of CR sources to a constant for R > iibk = 10 — 15 kpc. 



assume a halo height of 4 kpc and we set the CR source density profile to a uniform value 
beyond a given break radius. We found that 7-ray data point to a flat source distribution 
beyond R ^ 10 kpc. 

The solutions are not unique: alternative CR propagation models can be considered, 
with, e.g., a non uniform diffusion coefficient [10]. On the other hand, gni might be 
overestimated due to large amounts of dark gas in the outer disc of the Milky Way not 
accounted for by radio/mm- wave data [11]. 

* * * 

The Fermi LAT Collaboration acknowledges support from a number of agencies 
and institutes for both development and the operation of the LAT as well as scien- 
tific data analysis. These include NASA and DOE in the United States, CEA/Irfu and 
IN2P3/CNRS in France, ASI and INFN in Italy, MEXT, KEK, and JAXA in Japan, 
and the K. A. Wallenberg Foundation, the Swedish Research Council and the National 
Space Board in Sweden. Additional support from INAF in Italy and ONES in France 
for science analysis during the operations phase is also gratefully acknowledged. 



REFERENCES 

[I] DiGEL S. W. et al, Astrophys. J., 463 (1996) 609; 555 (2001) 12. 
[2] Atwood W. B. et at, Astrophys. J., 697 (2009) 1071. 

[3] Abdo A. A. et al, Astrophys. J., 710 (2010) 133. 

[4] ACKERMANN M. et ai, Astrophys. J., 726 (2011) in press; arXiv:1011.0816. 
[5] Abdo A. A. et ai, Astrophys. J., 703 (2009) 1249. 
[6] Strong A. W. et al, Astron. Astrophys., 537 (2004) 763. 
[7] Nakanishi H. and Sofue Y., Publ. Astron. Soc. Jpn., 58 (2006) 847. 
[8] Strong A. W., Moskalenko I. V. and Ptuskin V. S., Annu. Rev. Nucl. Part. S., 57 
(2007) 285. 

[9] Strong A. W., Moskalenko I. V. and Reimer O., Astrophys. J., 613 (2004) 962. 
[10] EvOLl C. et ai, J. Cosmol. Astropart. P., 10 (2008) 18. 

[II] Grenier la., Casandjian, J.-M. and Terrier R., Science, 307 (2005) 1292. 



